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Preface 

The  use  of  the  fluidized  solids  technique  was  developed 
very  largely  by  the  petroleum  and  chemical  industries,  for 
processes  where  the  very  high  heat  transfer  coefficients 
and  the  high  degree  of  uniformity  of  temperature  within  the 
bed  enabled  the  development  of  processes  which  would  otherwise 
be  impracticable. 

Fluidized  solids  are  now  used  quite  extensively  in  many 
industries  where  it  is  desirable  to  bring  about  intimate 
contact  between  solid  particles  and  gas  stream.  The  principles 
of  fluidization  are  applied  in  the  drying  processes  to  get  a 
f luidized-bed  drying  operation. 

The  purpose  of  this  study  is  to  develop  an  analytical 
model  of  the  f luidized-bed  drying  process  which  show  the 
behavior  of  the  capital  and  power  cost  of  such  a  system  in 
response  to  changes  in  fluidization  velocity. 

I  would  like  to  express  my  gratitude  to  my  thesis  advisor 
Professor  Robert  F.  Allen  for  his  most  valuable  advice  and 
guidance  during  this  study.  I  would  also  like  to  thank 
Professor  Milton  Franke,  my  reader,  for  his  help  during  the 
study. 

The  demanding  nature  of  any  thesis  project  is  hard  on  the 
student  accomplishing  the  thesis,  but  it  is  equally  as  hard  on 
the  spouse  of  the  student.  I  am  indebted  to  my  wife.  Dr. 
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MOUSSA  I.  M.  MOSTAFA 


ii 


Contents 


Page 

ii 


Preface . 

List  of  Figures . v 

List  of  Tables .  vii 

Nomenclature . ix 

I.  Introduction . „ .  1 

A.  The  Phenomenon  of  Fluidization .  1 

B.  Objectives . . . .’ .  4 

C .  Me  thodo  logy .  4 

II.  Engineering  and  Economic  Aspects  for 

Gluidized-Bed  Drying.... . 6 

A.  Engineering  Aspects .  6 

1 .  Fluidized-Beds . 6 

2.  Drying  Operation .  8 

3.  Compression  of  Air . 14 

4.  Heat  Exchange .  17 

B.  Economical  Aspects . 18 

III.  Economic  Model .  21 

A.  Continuous  Direct  Contact  Drying..... .  22 

1.  The  Dryer  Capital  Cost . . .  22 

2.  The  Air  Heater  Capital  Cost .  26 

3.  The  Steam  Cost .  27 

4.  The  Compressor  Capital  Cost .  28 

5.  Compressor  Power  Consumption  Cost .  32 

B.  Total  Cost . 32 

IV.  Presentation  and  Discussion  of  Model  Results .  38 

A.  Capital  and  Power  Cost...... . 42 

1.  Minimum  Fluidization  Velocity . . 

2.  Solid  Flow  Rate .  49 

3.  Other  Auxiliary  Factors . 53 

B.  Separate  Cost  Items .  56 


1.  Capital  Cost.... . 

2.  Power  Consumption  Cost 


56 

63 


Contents 

Page 


C.  Summary,  Conclusions,  and  Recommendations...  68 

Bibliography .  71 

Appendix  A:  Fibonacci  Algorithm..... . .  73 

Appendix  B:  FORTRAN  Code.. . 77 

Vita. . 81 


! 


ABSTRACT 

This  thesis  investigates  the  economic  characteristics 
of  a  f luidized-bed  drying  process.  It  focuses  on  the  way 
in  which  fluid  velocity  impacts  the  capital  and  power 
consumption  costs  of  such  a  system. 

The  engineering  and  economic  relationships  employed 
in  the  model  are  developed  and  a  sensitivity  of  the  optimal 
design  of  a  f luidized-bed  drying  process  based  on  capital 
and  operating  cost  is  explored  through  variation  in  the 
technical  parameters  of  the  model v 


viii 


NOMENCLATURE 


Hq  Moisture  content  of  leaving  air,  lb.  water/lb.  dry  air. 

i  Average  rate  of  return  on  invested  capital^ (frac./yr. ) 

im  Minimum  acceptable  return  rate  on  invested  capital. 

K  Proportionality  constant, 

k  Specific  heat  ratio  for  air,  dimensionless. 


ix 


L  Length  of  the  dryer,  ft. 

Lr  Supscript  for  large  size  equipment. 

M  Maintenance  charge, (frac . /yr. ) 

n  Constant, 

n^  Equipment  size-cost  exponent. 

TI2  Allowable  life  for  depreciation  of  equipment. 

P  Pressure,  lb^./sq.  in. 

P^  Entering  gas  pressure,  lb^/sq.  in. 

Pj  Compressed  gas  pressure,  lb^/sq.  in. 

0  Flow  rate,  cu.  ft. /min. 

Qa  Rate  of  heat  transfer  in  the  dryer,  B.t.u./hr. 

0H  Total  heat  transfared,  B.t.u./hr. 

r  Fixed  charge  factor/(frac./yr. ) 

R  Gas  constant,  ft./ (lb.  mole)  (°K) 

Re  Reynolds  number,  dimensionless. 

S  Supscript  for  small  size  equipment, 

t  Income  tax  rate, (frac . /yr . ) 

TDi  Inlet  air  dry  bulb  temperture,  °F. 

TDq  Outlet  air  dry  bulb  temperature,  °F. 

T  ^  Inlet  air  temperature,  °F. 

Tgo  Outlet  air  temperature,  °F. 

T  Steam  temperature,  °F. 

s 

Tw  Air  wet  bulb  temperature,  °F. 

umf  Minimum  fluidization  velocity,  ft. /sec. 


x 


ut  Terminal  velocity/  ft. /sec. 

U  Overall  volumetric  heat  transfer  coefficient/ 

B. t.u./ (hr. )  (cu.  ft.)  (of). 

U„  Overall  heat  transfer  coefficient/  B. t.u. /(hr.) 

M  (sq.  ft.)  (°F). 

w  Gas  weight  rate,  lb. /sec. 

W  Pounds  of  dry  solid,  lb. /hr. 

v  Gas  volume. 

V  Dryer  volume,  cu.  ft. 

x ^  Entering  solid  moisture  content,  lb.  water/lb.  dry  air. 

xq  Leaving  solid  moisture  content,  lb.  water/lb.  dry  air. 


Greek  Letter  Symbols 


AP  Pressure  drop,  lbf./sq.  in. 

AT„  Logarithmic  mean  temperature  difference  between 

air  and  steam,  °F. 


ATm  Logarithmic  mean  temperature  difference,  °F. 


y  Fluid  viscosity,  lb. /(ft.)  (sec.) 

x  Constant 

p  Air  density,  lb./cu.  ft. 

Pg  Gas  density,  lb./cu.  ft. 

p_  Particle  density,  lb./cu.  ft. 

<b  Particle  shape  factor,  dimensionless. 
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AN  INVESTIGATION  INTO  THE  RELATIONSHIP 
OF  FLUIDIZATION  VELOCITY  TO  THE  OPTIMAL  ECONOMIC 
DESIGN  OF  A  FLUIDI ZED-BED  DRYING  PROCESS 

I .  INTRODUCTION 

Fluidization  may  be  defined  as  the  phenomena  in  which 
the  gravitational  force  acting  orj  a  dense  swarm  of  particles 
is  counteracted  by  the  drag  of  an  upward  fluid  stream,  which 
cause  these  particles  to  be  kept  in  a  more  or  less  floating 
state  (1^)  . 

A,  The  Phenomenon  of  Fluidization 

The  term  fluidization  is  used  to  describe  a  certain  mode 
of  contacting  granular  of  finely  divided  solid  particles  with 
fluids.  The  physical  phenomena  occuring  in  the  passage  of 
fluids  through  unrestrained  particles  beds  is  described  by 
many  investigators.  Zenz  and  Othmer  (id)/  Leva  (11),  Kunii  and 
Levenspiel  (8) ,  Coulson  and  Richardson  (3) ,  Lewis  and  Bower- 
man  (13),  Ergun  (5)  and  Brownell  and  Katz  (2). 

Passing  a  fluid  upward  through  a  bed  of  fine  particles 
involves  certain  characteristic  states.  These  are  shown  in 
Figure  1  .  When  the  flow  rate  is  low,  fluid  percolates 
through  the  void  space  between  stationery  particles.  This  is 
the  state  of  a  "fixed  bed"  , 

•  With  an  increase  in  flow  rate,  an  "expanded  bed"  occurs. 
At  some  velocity  the  particles  are  all  just  suspended  in  the 
upward  floating  gas  or  liquid.  At  this  point  the  frictional 


force  between  a  particle  and  fluid  counterbalances  the  weight 
of  the  particle.  The  pressure  drop  through  any  section  of  the 
bed  approximately  equals  the  weight  of  fluid  and  particles  in 
that  section.  The  bed  is  referred  to  as  an  "incipiently 
fluidized  bed”  or  a  "bed  at  minimum  fluidization"  (8) . 

With  an  increase  in  flow  rate  beyond  minimum  fluidization, 
gas-solid  systems  behave  in  an  unstable  manner  with  bubbling 
and  channeling.  At  higher  flow  rates,  the  movement  of  solids 
becomes  more  vigorous.  Such  a  bed  is  called  an  "aggregative 
fluidized  bed". 

At  a  sufficiently  high  fluid  flow  rate  the  terminal 
velocity  of  the  solids  is  exceeded,  the  upper  surface  of 
the  bed  disappears,  entrainment  becomes  appreciable,  and 
solids  are  carried  out  of  the  bed  with  the  fluid  stream. 

The  fluid  bed  reactor  has  a  broad  usage  in  chemical 
engineering  industry  (8) .  However,  the  complex  nature  of  this 
fluid  bed  makes  accurate  mathematical  modeling  difficult. 

For  example.  Dale  (5)  describes  the  problem  of  waste  disposal 
in  the  Navy  in  the  late  sixties  and  the  Navy's  interest  in  the 
combustion  of  solid  material  in  a  f luidized-bed  reactor.  The 
lack  of  fundamental  design  knowledge  at  that  time  was  pro¬ 
hibiting  accurate  design  and  control  of  such  a  system.  Dale 
(4)  had  limited  success  in  building  a  mathematical  model  for 
non-catalytic  reactions  in  a  f luidized-bed  reactor. 

This  study  follows  the  spirit  of  the  Dale  research  in 

that  it  represents  a  further  exploration  of  an  analytical 

model  for  f luidized-bed  operations.  However,  the  goal  of  this 

3 


research  study  is  to  investigate  the  economic  characteristics 
of  a  f luidized-bed  drying  process.  In  particular,  this  study 
will  focus,  on  the  way  in  which  fluid  velocity  impacts  the 
capital  and  power  consumption  costs  of  such  a  system.  If 
successful  this  model  can  be  used  as  a  design  tool  in  studies 
of  f luidized-bed  drying  systems. 

B.  Objectives 

The  specific  objective  of  this  work  is  to  develop  an 
analytical  model  of  the  f luidized-bed  drying  process,  which 
will  highlight  the  behavior  of  the  capital  and  power  cost  of 
such  a  system  in  response  to  changes  in  fluidization  velocity. 
We  are  especially  interested  in  determining  the  extent  to 
which  minimum  system  cost  is  sensitive  to  variation  in  fluid¬ 
ization  velocity.  A  secondary  objective  of  this  study  is  to 
set  out  the  cost-velocity  relationships  of  each  of  the  major 
equipment  components  of  the  f luidized-bed  drying  system, 

C.  Methodology 

The  method  of  approach  used  for  this  investigation  is 
analytical.  Using  the  available  literature  on  chemical 
engineering  process  economics  and  cost  estimation,  a  model  is 
developed  relating  fluid  velocity  of  the  f luidized-bed  drying 
operation  to  the  capital  and  power  consumption  cost  of  the 
system. 

This  is  done  in  two  steps.  First  the  basic  engineering 
relationships  pertaining  to  heat  transfer,  pressure  change. 


equipment  volume  and  surface  area  of  equipment  are  specified 


so  as  to  highlight  the  role  of  fluid  velocity.  Second,  these 
engineering  parameters  are  related  to  capital  and  power 
consumption  costs  based  on  the  empirical  estimates  of  equip¬ 
ment  cost  size  relationships.  The  resulting  aggregate  cost 
function  represent  an  objective  function  to  be  minimized 
subject  to  a  velocity  range  constraint.  The  upper  and  lower 
bounds  of  this  velocity  constraint  are  the  terminal  and 
minimum  fluidization  velocities,  respectively. 

The  optimization  technique  used  to  solve  the  cost  model 
is  the  Fibonacci  search  method  I9) . 
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II.  ENGINEERING  AND  ECONOMIC  ASPECTS 


FOR  FLUIDIZED-BED  DRYING 


This  chapter  will  summarize  the  key  theoretical  engineer¬ 
ing  relationships  and  economic  cost  concepts  which  will  be  used 
in  developing  of  the  aggregate  cost  model  of  f luidized-bed 
drying  operations.  The  engineering  relationships  are  mainly 
analytical  equations  which  describe  the  physical  operation, 
while  the  primary  economic  cost  concepts  are  the  equipment 
cost-size  relationships  and  its  conversion  into  annual  costs 
by  use  of  an  annual  capital  charge  factor. 


A.  Engineering  Aspects 
Fluidized  Beds 

At  minimum  fluidization,  values  of  mass  velocity  and 
voidage  can  be  related  by  Leva  (3) . 


G,  =  °-°°5  DP  g  pq<ps  -  0g>  s  tk 

u'1  -  S.f> 


(1) 


Where  Gmf  =  fluid  superficial  mass  velocity  for  minimum 
fluidization,  lb. /(sec.)  (sq.  ff4  ; 

D  =  particle  diameter,  ft.; 

hr 

g  =  local  acceleration  due  to  gravity, 
32.2  f t . / (sec. ) (sec. ) 


Pg  =  gas  density,  lb/cu.  ft.; 


d  =  solids  density,  lh/cu.  ft>; 

5 

<t>s  =  particle  shape  factor,  dimensionless; 


em£  =  voidage  at  minimum  fluidization, 
dimensionless; 

y  =  fluid  viscosity,  lb. /(ft.)  (sec.) 
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If  and/or  4»g  are  unknown,  a  modification  made  by 

Wen  and  Yu  (15)  can  be  used  which  simplifies  the  expression 
of  minimum  fluidization  velocity.  Wen  and  Yu  found  that  for 
a  wide  variety  of  systems: 

— —  -  14  and  1  -emf  _11  (2) 

*semf  ,2  2 

.  *3  €mf 

which  when  substituted  in  the  minimum  fluidization  velocity 
expression  given  by  Kunii  and  Levenspiel  (6)  gives: 

umf  =  DP  (ps  ~  V  g  (3) 

1650  p 

where  u\ .  -  =  minimum  fluidization  velocity,  ft/sec. 

The  gas  flow  rate  through  a  fluidized  bed  is  limited  on 
the  one  hand  by  Umf  and  on  the  other  hand  by  entrainment  of 
solids  by  the  gas.  When  entrainment  occurs  these  solids  must 
be  recycled  or  replaced  by  fresh  material  to  maintain  stead- 
state  operations.  The  upper  limit  of\  the  gas  velocity  is 
approximated  by  the  terminal  or  free-fall  velocity  of  the 
particles  which  is  reported  by  Kunii  and  Levenspiel  (8)  to 
be  given  by: 

Ut  =  g  (ps  “  pa}  DP  for  Re<  0,4 
18y 

and 

u.  =  f4  (ps  “  pa)2J  1/3  Dp  for  0.4  <  Re  <  500  (5) 

t  X  225 

where  ufc  is  the  terminal  velocity,  ft/sec.; 

Re  **  DP  pg  Ut  is  the  Reynolds  number, 

y  dimensionless. 


Kunii  and  Levenspiel  (8)  calculated  the  ratio  of  terminal 


('•  ■■ 


velocity  ufc  to  minimum  velocity  um^  directly  from  the 

equations  of  minimum  fluidization  and  terminal  velocities. 

For  a  typical  fluidization  process  of  fine  solids  this  ratio 
is: 

=  91.6  (6) 

umf 

This  ratio  will  govern  the  range  of  velocities  that  will 
be  considered  feasible  in  the  analysis  of  the  economic 
cost  model  (17) ,  (12) ,  and  (ID) . 

Drying  Operation 

Drying,, as  defined  by  Perry  and  Chilton  (15),  generally 
refers  to  the  removal  of  a  liquid  from  a  solid  by  evaporation. 

In  their  excellent  survey  of  drying  operations  in  the 
chemical  industry.  Perry  and  Chilton  (15)  presented  a  class¬ 
ification  of  dryers  based  on  heat  transfer  which  is  shown  in 
Figure  2  . 

As  can  be  seen  in  Figure  2  ,  the  f luidized-bed  drying 
process  is  a  direct  contact  continuous  drying  operation.  This 
process  is  described  by  Coulson  and  Richardson  (3)  as  follows: 
"Heated  air,  or  hot  gas  from  a  burner,  is  passed  via  a  plenum 
chamber  and  a  diffuser  plate  into  the  f luidized-bed  of  material, 
from  which  it  passes  to  a  dust  separator.  The  diffuser  plate 
is  fitted  with  suitable  nozzles  to  prevent  back-flow  of 
solids.  Wet  material  is  fed  continuously  into  the  bed  through 
a  rotary  valve  and  mixes  immediately  with  the  dry  charge. 

Dry  material  overflows  via  a  downcomer  to  an  integral  after- 
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All  type*  of  dryers  used  for  producing  a  dry, 
♦olid  product  from  o  wot  food 


Diroct  dry  or* 

Hoot  tronsfer  for  drying  is 
accomplished  by  diroct  cent  act 
between  the  wet  solid  and  Hot 
gases.  The  vaporized  liquid  is 
carried  away  by  the  drying 
medium;  i.e,,  the  hot  gases. 
Direct  dryers  might  also  be 
termed  convection  dryers 


rodiont-heat  dryer*, 
dielectric-heat  dryer* 

The  operation  of  rodiant  heot 
dryers  depends  on  the  genera¬ 
tion,  transmission,  and  absorp¬ 
tion  of  mfrored  toys. 

Dielectric-heat  dryers  oper¬ 
ate  on  the  principle  of  heat 
generation  within  the  solid  by 
placing  the  latter  in  a  high- 
frequency  electric  field 


Indirect  dryer* 

Hoot  for  drying  i*  trans¬ 
ferred  to  the  wet  solid  through 
o  retaining  wall.  The  vaporized 
liquid  is  removed  independ¬ 
ently  of  the  heating  medium. 
Rate  of  drying  depends  on  the 
contacting  of  the  wet  moterial 
with  hot  surfoce*.  Indirect 
dryers  might  olso  be  termed 
conduction  or  contact  dryer* 


Operchon  is  continued  without 
interruption  os  long  as  wet  feed 
is  supplied.  It  is  apparent  that 
ony  continuous  dryer  can  be 
operated  intermittently  or  batch- 
wise  if  so  desired 


Dryers  oie  designed  to  oper¬ 
ate  on  o  definite  size  of  batch  of 
wet  feed  for  given  time  cycles. 
In  botch  dryets  the  conditions  of 
moisture  content  and  temperature 
continuously  change  at  any  point 
in  the  dryer 


Direct  continuous  type* 

1.  Continuous-tray  dryers  such  as 
continuous  metal  belts,  vibrating 
troys  utilizing  hot  gases,  verticoi 
turbodryers 

2.  Continuous  sheeting  dryers.  A 

continuums  sheet  of  material  posses 
through  the  dryer  either  os  festoons 
or  as  o  rout  sheet  stretched  on  a 
pin  frame 

3.  Pneumotic  conveying  dryer*. 

In  this  type,  dryng  is  often  done 
in  coniunction  with  grinding.  Mo* 
tenal  conveyed  in  htgh-tempero- 
tire  high- velocity  gases  to  e  cy¬ 
clone  collector 

4.  Rotary  dryers.*  Material  is  con¬ 
veyed  and  showered  inside  o  ro¬ 
tating  cylinder  through  which  hot 
gases  flow 

5.  Spray  dryers.  Dryer  feed  must 
d e  coooble  of  otomization  by 
either'  a  centrifugal  disk  or  o 
nozzle 

6.  T h re ugh-circulotien dryer*.  Mo- 
tenal  is  held  on  o  continuous  con¬ 
veying  screen,  and  hot  oir  is 
blown  tmough  it 

7  Tunnel  dryers.  Material  on  trucVs 
is  moved  through  o  tunnel  m  con¬ 
tact  w.th  hot  gases. 

8.  fluid  beds.  Solids  ore  fluidized 
m  a  stationary  tom,.  May  also 
hove  indirect-heat  coils 


Continuous 

Drying  is  accomplished  by 
moterial  passing  through  the 
dryer  continuously  and  in  con¬ 
tact  with  o  hot  surioce 


Botch  indirect  dryers  are  gen¬ 
erally  well  adapted  to  operate 
under  vocuum.  They  may  be  di¬ 
vided  into  agitated  and  non- 
ogitated  types 


Direct  batch  typac 

1.  Cylinder  dryer*  for  continuous 

1.  Batch  through-circulation  dry- 

sheets  such  as  pope',  eeilophone , 

ers.  Material  held  on  screen  bat- 

textile  piece  goods.  Cylinders  ore 

tom  troys  through  which  hot  air 

generally  steam-heated,  ond  ro- 

is  blown 

fate 

2.  Tray  ond  comportment  dryer*. 

2.  Drum  dryers.  These  may  be 

Material  supported  on  troys  which 

heated  by  steom  or  hot  water 

may  or  may  nor  be  on  removable 

3.  !  crew- convey  or  dryer*.  Al-  j 

trucks.  Air  blown  across  material 

though  th»se  dryers  ore  contm- 

on  frays 

ucus,  operation  under  a  vocuum 

3.  Fluid  bed*.  Solids  are  fiuidizod 

is  feo»ib<e.  Solvent  recovery  with 

in  a  stationary  cart  with  dust  fll- 

drying  is  possible 

ter  mounted  obove 

4.  Steam-tube  rotory  dryer*. 

Steam  ar  hot  wa»er  ten  be  used. 

Op*  rot.cn  on  slight  negative  pres- 

•ore  is  fsosibie  to  permit  solvent 

recovery  with  drying  if  des.red 

S.  Vibroting-tray  dryer*.  Heating 

accomplished  by  steam  or  hot 

« 

wo*er 

A.  Spuciol  type*  such  os  a  contm- 

unus  fabric  belt  moving  in  ctose 

contoc*  with  a  st  tom- heated 

platen.  Mote^ol  to  be  dried  lies 

on  the  bzlt  ond  receives  heat  by 

•  Certain  rotory  dryer:  may  be  o 

contact 

combination  of  indn*ct  end  drect  types; 

1.  Agitated-pan  dryers.  These  may 
operate  atmospherically  or  under 
vocuum,  and  con  handle  smaR 
production  of  nearly  any  form  of 
wet  solid,  i.e.,  liquids,  slurries, 
pastes,  or  granular  solids 

2.  Freeze  dryers.  Material  is  frozen 
prior  to  drying.  Drying  in  frozen 
state  is  then  done  under  very  high 
vacuum 

3.  Vocuum  rotary  dryer*.  Mate¬ 
rial  is  ogitated  in  o  horizontal, 
str.hcnery  shell.  Vocuum  may  not 
olways  be  necessary.  Agitator  may 
be  steam-heated  in  addition  to  the 
shell 

4.  Vocuum-troy  dryer*.  Heating 
done  by  contact  with  steam- 
heated  or  hot- water -heated  shelves 
on  which  die  material  Net.  No 
ogitotion  involved 


e.g.,  ho*  gases  first  heo*  on  inner  shell 
and  then  pass  between  an  inner  and 
outer  shell  in  contact  with  the  wet  solid. 


a*0,  v, gV 
»*•  '-5 


»v  ,e? 


Figure  2  Classification  of  dryers  ba$ed  on  heat  transfer. 
Perry  and  Chilton  (15) 
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cooler".  Figure  3  shows  a  schematic  flow  sheet  for  a 
typical  equipment  configuration  for  a  fluidized  bed  drying 
opertion. 

Since  the  f luidized-bed  drying  operation  is  mainly  the 
removal  of  moisture  content  from  solid  particles  by  direct 
contacting  with  hot  air,  it  requires  the  basic  equipments 
for  heating,  blowing  of  air  and  fluidizing  of  solid  particles. 

As  shown  in  Figure  3  ,  air  enters  the  heat  exchanger 
with  the  following  specifications: 
moisture  content, 

dry-bulb  temperature,  TDi 

wet-bulb  temperature,  T^ 

Pressure,  P^ 

The  steam  enters  the  heat  exchanger  with  a  temperature 
Tg/  and  increases  the  air  temperature  to  Tg0* 

The  function  of  the  compressor  is  to  increase  the  air 

pressure  from  P^  to  P2  to  obtain  the  pressure  drop  required 

for  fluidization  of  solid  particles.  In  the  fluidization 

column,  the  direct  contact  of  heated  air  with  the  wet  solid 

results  in  the  removal  of  solid  moisture  content  from 

initial  quantity  X^  to  final  quantity  X^.  As  a  consequence 

the  moisture  content  of  the  air  is  increased  from  H.  to  H  . 

i  o 

The  design  used  in  this  research  of  the  fluidized  bed 
drying  process  is  an  approximation  of  the  direct  contact 
rotary  dryer.  The  heat  transfer  mechanism  employed  in  such  a 
system  is  given  by  Perry  and  Chilton  (15)  as: 


10 


where  Qa  =  total  heat  transferred,  B.t.u/hr.; 

U  =  volumetric  heat  -  transfer  coefficient, 

B.t.u./(hr)  (cu.  ft  dryer  volume)  (°FJ ; 

V  “  dryer  volume,  cu.  ft.; 

AT  =  true  mean  temperature  difference  between  the 
m 

hot  gases  and  material,  °F. 

When  a  considerable  quantity  of  surface  moisture  is  removed 
from  the  solids  and  the  solids  temperatures  are  unknown,  a 
good  approximation  of  (At^)  is  the  logarithmic  mean  between 
the  wet-bulb  depression  of  the  drying  air  at  the  inlet  and 
exit  of  the  dryer. 

Data  for  evaluating  U  were  shown  in  Perry  and  Chilton 
(15)  and  take  the  form 

U  =  K  Gn/D  (8) 

a 

where  K  =  a  proportionality  constant; 

G  =  gas  mass  velocity,  lb/ (hr.) (sq.ft.)  dryer  cross- 
sectional  area. 

D  =  dryer  diameter,  ft >; 
n  =  a  constant 

McCormick,  as  reported  by  Perry  and  Chilton  (14)  compared 
all  available  data  and  concluded  the  value  of  n  =  0.67  is 
probably  the  most  truly  representative  of  commercial  equipment 

The  drying  process  is  related  to  psychrometry  which  is 
concerned  with  determination  of  the  properties  of  gas-vapor 
mixture.  The  air-water  vapor  system  is  by  far  the  one  most 


CP 


commonly  encountered.  Terminology  and  relationships  pertinent 

to  psychrometry  as  reported  by  Perry  and  Chilton  (15)  are: 

Absolute  humidity  H  equals  the  pounds  of  water  vapor  carried 

by  one  pound  of  dry  air. 

Percentage  absolute  humidity  is  defined  as  the  ratio  of 

absolute  humidity  to  saturation 
humidity. 

Percentage  relative  humidity  is  defined  as  the  partial  pres¬ 
sure  of  water  vapor  in  air  divided  by 
the  vapor  pressure  of  water  at  the  given 
temperature. 

Dew  point  or  saturation  temperature  is  the  temperature  at 

which  a  given  mixture  of  water  vapor 
and  air  is  saturated 

Humid  heat,  Cg  is  the  heat  capacity  of  one  pound  of  dry  air 


and  the  moisture  it  contains. 


C  =  0.24  +  0.45  H 
s 


■:  v-  v 


Wet-bulb  temperature:  is  the  dynamic  equilibrium  temperature 

attained  by  a  water  surface  when  the 
rate  of  heat  transfer  to  the  surface 
by  convection  equals  the  rate  of  mass 
transfer  away  from  the  surface. 

Psychrometric  charts  for  air-water  vapor  systems  have 
been  developed  and  are  used  to  find  the  properties  of  a  system 
at  given  dry-bulb  and  wet-bulb  temperature. 

Equation.  (7)  will  be  used  for  calculation  of  the  dryer 
volume  and  equation  (8)  will  be  used  for  determination  of 
dryer  diameter. 
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Compression  of  Air 

Pressure  drop  is  required  for  the  flow  of  air  through 
the  beds  of  granular  solids.  It  varies  due  to  the  variety  of 
granular  materials  and  their  packing  arrangement.  This 
pressure  drop  can  be  created  in  the  f luidized-bed  drying  system 
by  either  compressors  or  fans.  Fans  are  used  for  low  pres¬ 
sure,  in  general,  for  pressure  drop  of  less  than  0.5  lb^/sq.  in 
(15) .  Centrifugal  compressors  are  widely  used  to  handle  large 
volumes  of  gas  at  pressure  rises  from  0.5  up  to  several 
hundred  lb^/sq.  in.  The  centrifugal  compressor  cost  estimation 
by  Happel  and  Jordon  (7)  is  given  by  equation  (48)  and  for 
fans  are 

Purchased  cost  =  6.7  (cfm)  0,88  $ 

Installed  cost  =  2.78  (6.8)  (cfm)  0,68  f 

Where  cfm  is  the  cubic  feet  per  minute  gas  flow  rate 
comparing  this  fan  cost  estimation  equations  with  the 
centrifugal  compressor  cost  estimation  in  equation  (48) 
reveals  that  compressors  will  be  more  expensive  than  fans. 

Due  to  the  limitation  of  fans  in  a  fine  particles  fluidized- 
bed  this  study  will  use  the  centrifugal  compressor. 

In  any  continuous  compression  process  the  relation  of 
absolute  pressure  P  to  volume  V  is  expressed  by  the  formula 

PVn  =  K  (9) 

The  plot  of  pressure  vs.  volume  for  each  of  exponent  n  is 
known  as  a  polytropic  curve.  Since  the  work  Wp  performed  in 
proceeding  from  P^  to  Pj  along  any  polytropic  curve  as  in 
Figure  4  is: 


it  follows  that  the  amount  of  work  required  is  dependent  upon 
the  polytropic  curve  involved  and  increases  with  increasing 
values  of  n.  The  path  requiring  the  least  amount  of  input 
work  is  n  =  1,  which  is  equivalent  to  isothermal  compression. 
For  adiabatic  compression,  i.e.  no  heat  added  or  taken  away 
during  the  process,  n  =  k  =  ratio  of  specific  heat  at  constant 
pressure  to  that  at  constant  volume. 

The  adiabatic  head  is  calculated  by: 

k  k-1  _ 

tt  _  **  T»*m  r  /  n  \  V  ■  1  1  /  1  1  \ 


ad  k 


5l  R  T!  C  <p2.) 


Where  R  =  gas  constant  =  1544  f  t./ ( lb.  mole)  (°iy  ; 

T  =  absolute  temperature,  °R. 

Since  the  work  performed  on  the  gas  is  equal  to  the  product 
of  the  head  and  weight  of  gas  handled,  the  adiabatic  horsepower 


HP  = 


k  WRTi  r,  P2.  k  -  11 

“k-Y  • -  I'pi'  1 

K  1  550  ± 


k-1 

HP  =  0.00436  0  P.  (  k  )  f(|y)  k  -1] 
'  1  k-1  L 


where  w  =  weight  of  gas,  lb/sec;? 

Ha(j  =  adiabatic  head,  ft>; 

0  =  flow  rate  cu.  ft/rain.; 

P^  =  pressure,  lbf./sq.ia 

Equation  (13)  will  be  used  for  compressor  capital  cost  as 
well  as  for  power  consumption  cost  calculations  in  developing 


cost  model 


Heat  Exchange 

The  basic  design  equation  for  the  heat  exchanger  stated 
by  Perry  and  Chilton  (15)  is: 


(14) 


Where  dAH  is  the  element  of  surface  area  required  to  transfer 
an  amount  of  heat  dQ  at  a  point  in  the  exchanger  where  the 
over-all  bulk  temperature  difference  between  the  two  streams 
is  Th. 


Equation  14  can  be  integrated  to  give  the  outside  area 
required  to  transfer  the  total  heat  load  Q 

H 


(15) 


The  heat  exchanger  used  in  the  design  of  drying  process 
uses  the  steam  to  heat  the  entering  air.  If  we  neglect  the 
heat  losses  to  the  surroundings,  an  energy  balance  about 
the  heat  exchanger  gives 


,heat  gained.  _  ,heat  lost, 

by  air  ‘  “  lby  steam  '  (16) 

For  conditions  normally  met  in  practice  these  continuous 
operations  can  be  considered  equilibrium  processes,  with  exit¬ 
ing  gas,  and  exiting  steam  all  at  the  same  temperature.  Thus 
with  F  (lb.  dry  gas/hr)  as  the  feed  rate  of  air,  the  energy 
balance  becomes 


F  Ce  (T  -  T  . )  =  U  A  AT„ 


(17) 


Where 


Tg  is  the  temperature  of  the  steam 
T  ^  is  the  temperature  of  the  air 
Cg  is  the  humid  heat  capacity  of  1  lb.  of 
dry  air  and  the  moisture  it  contains. 

For  most  engineering  calculations# 

C  =  0.24  +  0.45H 
s 

Where  0.24  and  0.45  are  the  heat  capacities  of  dry 
air  and  water  vapor,  respectively  and  both  are  assumed 
constant  (15) . 

Equation  (14)  is  used  in  estimating  the  capital  cost 
of  the  heat  exchanger  and  equation  (17)  will  be  used  in 
determining  the  cost  of  the  steam  in  developing  the 
cost  model. 

B.  Economical  Aspects 

The  principal  economic  aim  of  any  company  as  stated 
by  Freemantle  16)  is  the  efficient  use  of  available 
resources  in  order  to  maximize  benefit  over  the  agreed 
time  span.  Happel  and  Jordan  (7)  discussed  in  detail  the 
chemical  process  economics  which  provide  a  working  tool 
to  assist  in  applying  technical  information  to  the  economic 
design  and  operation  of  chemical  process  plants. 

The  essential  economic  relationship  of  interest  is 
the  equipment  cost  -  size  relationships  which  relates 
the  purchased  cost  of  equipment  to  the  size,  volume,  or 
power  of  the  equipment. 

The  general  form  of  the  equipment  cost  -  size 


relationship  is 
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Cost  =  k  (A) n 1 


(18) 


or 


(Cost)  Lr  /^Lr\nl 
(Cost)  s  -  [A sJ 


(19) 


where  A  =  some  characteristic  size  measurement 

such  as  volume,  area,  or  horsepower, 
k  =  a  constant,  the  value  of  the  cost 

when  A  is  unity. 

Lr  and S  =  subscripts  designating  larger  and 
smaller  sizes. 

Extensive  studies  of  process  equipment  exist  which 
generally  report  a  finding  of  0.6  or  0.7  for  the  value  of 
nl  (1),  The  evidence  is  so  favorable  to  this  value  that 
the  so  called  "point  6  value"  has  emerged  among  process 
engineers  to  describe  the  cost  -  size  relationship.  That 
is,  in  general  one  expects  a  10  percent  increase  in  equipment 
size  to  yield  a  6  per  cent  increase  in  capital  costs.  The 
economic  impact  of  an  nl  value  less  than  unity  is  that  there 
are  economy  of  large  size  operation  available  for  many 
engineering  processes. 

The  equipment  capital  cost  can  be  transformed  to  an 
annual  capital  -  related  fixed  charge  by  use  of  a  fixed 
charge  factor  r  which  is  defined  by  Happel  and  Jordan 
(7)  as: 


1 


i  -  0.35  1 


where 


i  =  minimum  acceptable  return  rate  on 
m 

invested  capital,  (frao/yr.)  ; 
n 2  =  allowable  life  for  depreciation  of 

equipment; 

=  straight-line  depreciation  factor; 
i  =  average  rate  of  return  on  invested 
capital,  (fracVyrJ  ; 

t  =  income  tax  rate,  Federal  plus  state, 
(frac/yrj  ; 

m  =  maintenance  charge,  (fracyyrO  ; 

While  the  fixed  factor  charge  will  vary  equipment 
and  firms,  Happel  and  Jordan  argue  that  a  value  of 
r  =  0.5  is  a. reasonable  approximation  for  most  engineering 
applications. 


III.  ECONOMIC  MODEL 


There  are  two  fundamental  economic  engineering 
problems.  One  consists  in  simply  finding  the  optimum 
capital  cost  for  the  equipment  to  accomplish  the  desired 
function.  Such  a  situation  exists  for  example  in  choosing 
the  optimum  diameter  for  tanks  and  vessels.  The  other 
consists  in  balancing  incremental  investment  costs  against 
possible  operating  savings.  In  fluid  flow  applications 
for  example,  one  might  balance  the  cost  of  larger  diameter 
pipe  against  the  saving  in  pumping  costs,  or  in  problems 
in  heat  transfer  one  may  balance  the  cost  of  extra  ex¬ 
changer  surface  against  the  cost  of  the  possible  heat 
saving. 

The  research  work  focuses  on  the  optimum  capital 
cost  for  the  equipment  of  a  continuous  direct  contact 
drying  process  as  one  of  the  physical  applications  of  the 
fluidized  bed  processes. 

In  studying  the  economics  of  process  equipment  one 
usually  assumes  a  constant  value  for  all  the  engineering 
variables  involved  except  one.  The  total  cost  function 
is  then  set  up  in  terms  of  this  variable  and  employs  one 
of  the  optimization  techniques  to  find  an  optimum  value. 
The  optimization  technique  that  will  be  used  in  this  study 
is  the  Fibonacci  search  method  (9)., 

The  main  feature  of  this  model  is  that  according  to 
the  calculation  of  the  dryer  heat  load  from  the  solids 
feed  rate  and  the  moisure  to  be  evaporated,  the  gas  flow 


rate  which  is  a  function  of  the  gas  velocity  is  calculated. 
From  the  known  dryer  heat  load,  head  -  transfer  coefficient 
and  the  logarithmic  mean  temperature  difference,  the  volume 
of  the  dryer  can  be  calculated.  Relating  the  crossectional 
area  of  the  dryer  and  the  height  to  the  gas  velocity  will 
develop  the  economics  of  the  model  in  terms  of  gas  velocity 

A.  Continuous  Direct  Contact  Drying 

As  discussed  in  Chapter  II,  drying  is  accomplished  by 
the  transfer  of  sensible  heat  from  a  gas  to  a  solid. 

The  evaporation  of  liquid  from  the  solid  and  removal 
of  the  evaporated  liquid  as  a  vapor  in  the  exit  gas.  Thus, 
a  dryer  is  both  a  heat  -  transfer  and  a  mass  -  transfer 
device  (8) . 

The  total  cost  of  the  drying  operation  consists  of  the 
capital  cost  of  the  dryer,  the  air  heater,  and  compressor 
and  the  operating  cost  of  steam  and  electric  power 
consumption.  We  will  discuss  each  separately  and  then  sum 
all  these  items  to  get  the  total  cost  model. 

1.  The  Dryer  Capital  Cost 

An  overall  humidity  or  water  balance  gives: 

W  (%i  -  =  F  (Hq-  H±)  (21) 

where  W  =  pounds  of  dry  solids  to  be  treated  per  hour; 

XQ  =  leaving  solid  moisure  content  ( lb,  water/lb.  dry 
solid) ; 


X^  =  entering  solid  moisure  content  (lb,  water/ lb. dry 
solid) ; 


Hq  =  humidity/  or  moistuje  content,  of  leaving 
gas  (lb.  water/ lb.  dry  gas); 

humidity,  or  moistuie  content,  of  entering  gas 
(lb.  water/lb.  dry  gas) ; 

F  *  m  dry  gas  /  hr. 


Thus, 


F  = 


W(xt  -  X0> 
(H0  -  H±) 


But  since  F  (lb  dry  gas/hr)  =  Ac  (cross  sectional  area  sq.  ft. 


Then 


x  p  (density  of  gas  lb.  )  x  u  (velocity  ft) 
y y  y  \  ca.  ft .  hr. 

F  =  “  Xoi  =  A  o  u 


(X0  -  H±) 


(22) 


Since  G,  the  mass  rate  of  flow  of  the  gas  in  the 
dryer  lb/ (hr.)  (sq.ft)  =  p  u 
from  Equation  (22)  we  have 


i,e. 


F  =  A_  0  u  =  A  G 
c  c 


_  F 


But  A  (cross  sectional  area  =  ?  u 
c  4 

Where  D  is  the  diameter  of  the  dryer 

‘4  F  -.0.5 


Thus 


or 


» -If  §  3' 


r 4  »«!  -  V  -l0-5 
D  "U-  p  W  (Hc  -  Hi)  J 


(23) 

(24) 


The  rate  of  heat  transfer  in  the  dryer  is  given  by: 


Oa  =  Ua  V  ATm  (25) 

Where  Qa  =  the  rate  of  heat  transfer  in  the  dryer  (B£ju/hr) 

U  =  the  overall  volumetric  heat  -  transfer  coefficient 
Si 


» 
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Equation  (25) 


Substituting  f or  Qa  f  rom  Equation  (27)  and  for  U 

C 

from  Equation  (29)  we  get 

1000  W  (X.  -  X  )  F0*5 

r,  _  '  1  _ O 


17.725  (p  u) 1,26  AT 

n 

56.417  W  (X±  -  XQ)  FC 
(p  U)1*26  AT 


Since 


W(X1  -  V 

Ho  -  Hi 


then 


rw(xi  •  V 

56.417  W  (X.  -  Xo)[  -H— J"Hg 
_ o _ i _ 

(p  u)1-26  4^ 


56.417  W(X±  -  XQ) 
1*0  -  H^0-5  (p  u  ) 


,1.26 


The  installed  cost  of  the  dryer  may  be  expressed  as 


a  function  of  its  peripheral  area,  A^  =  tt  D  L , 


Since 


where  L  is  the  length  of  the  dryer  expressed  in  feet. 

2  V  =  -7  D2  L 

4 


A  =  ir  D  4  V 

P  - 2 

n  D 


Substituting  the  expressions  for  V  and  D  from 


Equations  (24)  and  (31)  in  the  above  equation  for  A  gives 
56.417  [w(X±  -  XQ)  J1*5  ptt  p  u  (H  -  H.)  "1° ' 


P  {Ho  “  Hi)0*5  (  P  u  )1#26  AT 


4  W(X±  -  Xo) 
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From  the  correlations  of  the  costs  of  dryers  as 
presented  by  Happel  and  Jordan  (15)  ,  it  can  be  found  that 
the  installed  cost  of  a  carbon  steel  dryer  is 
Purchased  cost  =  195  (A  $ 

Installed  cost  =  3.0  (195)  $  (3 

P 


Where  A  is  the  peripheral  area  of  the  dryer 

Sr 


The  installed  cost  of  the  dryer  is  then 

^200.0  W  (Xi  -  XQ)  -|0,8 
l.  .  0.76  J 


Cost  =  585 


(3 


(p  u)  °*76  AT 


m 


The  annual  fixed  charges  on  the  dryer  is 

0.8 


Cost  =  r  ^  585  1 


l 


(p 


r  200  w  (x±  -  xQ)  y'°  \ 
u)0-76  AT  1  J 


(3 


m 


where  r  is  annual  fixed  charge  factor 


2.  The  Air  Heater  Capital  Cost 

The  heater  required  for  the  heating  of  the  inlet  air 
may  be  assumed  to  be  a  standard  steam  -  heated  exchanger 
of  the  usual  shell  and  tube  type,  and  made  of  carbon  steel 
From  the  energy  balance  of  equation  (21)  the  heating 
load  will  be 


W  (X  —  x  )  r 

Oh”  h*"-  h?  I0'24  (Tgo-  Tgl>  °-45  Bi(V  VJ 

Tgo  =  outlet  gas  temperature  (°F,) 

T  .  =  inlet  gas  temperature  (°F.) 


(3 


where 


°H  UH  AH  ATH 


(38) 


t? 


where 


U 


H  =  the  overall  heat  -  transfer  coefficient 


BiW  (hi;)  (sq.  f  tj  (°EJ 


ATh  = 


(T  -  T  .  )  -  (T  -  T  ) 


lai 


321 


Ln 


,T  -  T  . 
(  s  gi 


F.,  the  mean 


T  -  T 
x  go 


logarithm  temperature  difference  between  the  air 
and  steam 

Apj  =  the  peripheral  area  of  the  heater  (sq.  ftj  from 
Equations  (37)  and  (38) 


W(X.  -  xo) 

H  (H  -  H.)  Ut.  AT 


[0.24 (T 


-  T  .)  + 
go  gi 


i'  K  “‘H 

°-45  Hi  <Tgo  -  T,i0 
The  purchased  and  installed  costs  of  heat  exchangers 


(39 


have  been  estimated  by  Happel  and  Jordan  (15)  ass 

0.62 


Purchased  cost  =  0.5  (A„) 

rl 


Installed  cost  =  346  (A„) 

rl 


0.62 


$ 

$ 


Thus  the  annual  fixed  charge  on  the  air  heater  will  be 

Tir.  f  W(Xi  -  Xq)  [  0 . 24  (T  -  T  .) 

346  ' (H0  "  Hi)  UH  A  th  1  9  9 

,\  1 

$/yr.  (40 


0.45  H,  (Tgo  -  Tgi) 


3.  The  Steam  Cost 

The  cost  of  the  steam  used  to  heat  the  air  before  it 

enters  the  dryer,  assuming  24  hr.  x  330  day  =  8000  hours 

day  year 

of  operation/year  will  be 
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r  W(X .  -  X  )  -i 

Cost  =  8,000  Cs  [  _  H°  J  [0.24  (Tgo  -  Tgl)  + 

0.5  H,  (Tgo-Tgi)J  <4 

Where  C  =  cost  of  steam  $/10^  BJtji.  Happel  and 
s 

Jordan  (7)  estimated  C  to  be  0.8  $/10^  BAiu  in  1975 

s 

and  we  will  assume  the  energy  cost  is  increased  three 
times . 

The  gases  leaving  the  dryer  will  be  at  a  fairly  low 
temperature  and  will  be  assumed  to  have  no  value. 

4.  The  Compressor  Capital  Cost 

Power  consumption  is  an  important  factor  in  any 
process  using  fluidized  beds.  Figure  (5)  shows  the  power 
required  to  supply  sufficient  incoming  fluidizing  gas  at 
pressure  P3  is  that  necessary  to  compress  that  quantity  of 
gas  from  pressure  P3  to  Pj. 


Figure  5 .  Pressure  at  various  locations  in  fluidized  bed. 
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This  is  found  by  taking  a  mechanical  energy  balance  across 
the  compressor.  Thus  the  shaft  work  required  per  mole  of 
gas  compressed,  as  stated  by  Kunii  and  Levenspiel  (&)  is: 

,P 2  .  u?  -  u 

-  W_  «  ~  AL  +  /  •”  +  — ~ —  +  (friction  loss)  ( 

s  gc  •'p  p  2  gc 

For  adiabatic  reversible  operations  with  negligible 
kinetic  and  potential  energy  effects,  this  reduces  to  the 
ideal  work. 

The  adiabatic  head  is  calculated  as 

“ad  "  &1  RT1  [  <pf)  k  -  1  1  < 


where 


k  =  specific  heat  ration  =  1.4 


R  =  gas  constant  (1544  f  tV  ( lb.  mole)  (  RJ 
T  =  absolute  temperature  ,  °R. 

Since  the  work  performed  on  the  gas  is  equal  to 
product  of  the  head  and  the  weight  of  gas  handled,  the 


adiabatic  horsepower  is:  (11) 


HP  = 


‘  ad 
550 


Ji_  111  1[(!2}  TT  .  x  ] 

k-1  550  L1P,J  1  J 


HP  =  0.2616  0  P1  (j^)  [  (^) 

where  W  =  weight  of  gas  lb*/sec. 


2,  k 


=  adiabatic  head,  ft.; 

0  =  flow  rate,  cu.ftv'hr.; 


P  =  pressure  lb^./sq.  in. 

The  flow  rate  0  is  related  to  the  velocity  u  by 


the  relation 


AP  =  P3  -P4  «  P2  -  Px 

but  AP  is  calculated  for  the  fluidized  bed  ass 

AP  =  y44  ( 1-e)  (ps  -  pg)  J  (lbf./sq.  in.)  (50) 
where  L  =  bed  height,  ft./ 


p  =  gas  density,  lb./cu.ft.? 

p  =  solid  density,  lb./cu.ft.; 
s 

e  =  bed  voidage  fraction,  dimensionless/ 

g  =  acceleration  of  gravity,  32.2  f t . / (sec . ) (sec . ) ; 

gc  =  conversion  factor,  lb.  ft . / ( lbf . ) (sec . ) (sec . ) . 

P^  is  usually  the  atmospheric  pressure  =  14.7  lb^./sq.  in. 
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Thus 


AP  = 


U  -  e>  (Cs  -  Cg)a=f  P2  -  Pj 


144  (1  *  £)  <Ps  '  Pg*2“  P2  *  14'7 

gc 

P2  =  14*7  +  144  (1  "  e)  (ps  “  Pg)2 
4  V 

Since  L  =  - r,  substituting  for  V  and  D  from 

n 

Equations  (24)  and  (31) /  we  get 

4  (56.417  )  W  (X:  -  XQ)  F0,5 

L  "  *  Tt  »  )°'66  4T  (it-) 
v  >  m  \  tt  p  u  / 


0.34 


56.417  W  (X±  -  XQ)  ( p  u  )  °* 

W  (X.  -  X  ) 

AT  -  . . -1 . .  °.. 

m  <Ho  -  Hi) 

56.417  (0  U)  0,34  (H  -  H.) 
_  _  o  i 


Equation  (51)  becomes 
P2  =  14.7  +  0.4  • 


(P  u)0,34  (1  -  £)  (Pe  -  P)  g 

_  _  _  _  S 


Substitute  for  P~  in  Equation  (29) 

“  t  9.1  / 


Installed  cost  =  5873.68 


W  (X.  -  X  )  (£  ,) 

l  o  k-1 

p  (H  -  H.) 

K  O  1 


'  (p 2  expression)^  l)0,8 

L  14.7  “  1JJ  * 

The  annual  fixed  cost  of  the  compressor  is 

,  iazJ 

\W(X.  -  X  )  .  r  (P_  expression)  k 

r  =  5873 . 68|  -  ^  (i,)  l  -^477 - 

-TO. 8 

j  $/yr. 


31 


5.  Compressor  Power  Consumption  Cost 


The  cost  of  power  consumption  is  calculated  as 
follows: 


1.  The  annual  hours  =  8,000  hrk 

2.  The  horsepower  of  the  compressor  =  HP  from 
Equation  (57) 

3.  Efficiency  is  taken  to  be  0.75  (8) 

4.  The  cost  of  electrical  power  Ce  is  taken  as 
$0 . 24/kw.  hr.  since  it  was  estimated  in  1975  by 
Happel  and  Jordon  (7)  as  $0.08/Kw/hr. 

The  power  consumption  cost  =  8,000  (HP.)  (Ce)  $ 


=  41018.88  Ce 


0.75 


k-1 


W(Xi  -  xo) 

P  <HC  -  Hi) 


(j“i)  r^2  exPressionJ  k 


G 


14.7 


-1 


B.  Total  Cost 

The  total  capital  cost  is  obtained  by  adding  the  cost 
of  each  equipment,  represented  by  Equations  (36) ,  (40) 


and  (55) . 


J  r  200  W  (X.  -  X  K 

Capital  Cost  =  r  <585  - TTWc. -  ] 

l  1 (p  u)0*76  A  T_ 


0.8 


m 


annual  fixed  charge  on  the  dryer 


+  r 


{346  ( 


w<xl  -  V  [O-24  <Tgo  -  V  +  0.45  (Tgo 


T  .  )T 


(Ho  -  Hi)  uh  ath 


annual  fixed  charge  on  the  heater 


+  r  1  5873.68) 


W(X.  -  X) 


x  t( 

(k-l)l 


expression, 


k-l 


(k-l)i-  14,7 


0.8  . 

-•))  | 


annual  fixed  cost  of  compressor 


(57) 


The  total  operating  cost  is  considered  as  the  cost 


of  the  steam  used  to  heat  the  air  before  entering  the 
dryer,  and  the  electrical  power  consumption  by- the  compressor 


Operating  Cost  =  8,000  C 


W  (X,  -  Xj 


s  L  (Ho  -  H.) 


pi  0.24  (T  -T. )  + 
—  go  gi 


0.45  H.  (T  -  T  . ) 1 
i  go  gi  J 


go  gi 
+  41018.88  Ce 


W  (X^Xo>  (A) 


p  (Hq  -  Hi)  'k-l; 


,  k- 1 

r  (P2  expression)  -i 

*■  1477  "  1-’ 


(58) 


Thus  the  total  cost  =  Capital  cost  +  Operating  cost  or 

200  W  (Xi-XQ)  'l0*8 


Total  Cost  =  r 


{S85  ( 


(p  u)0,76  AT 


r] 


m 


+  r 


(W  [°-24  (Tgo'V  + 
V  gH  ath 


.0.62  , 

0.5  (Tgo-Tgl)]  J  j 
w  (X.-X J 


+  r 


5873 


(W 
.68  l- 

\  D 


1  ~o;  (r^r) 


p  (H0-Hi)  'k-l' 
k-1  0.8 


[(P2  expression)  —  .  J  j 

/  W  <XH-X  )  rn 

- 8'000  csl— -H-rr  1 

O  i 


24  (T  -T  . )  + 

go  gi 


°*45  Hi  (Tgo  -  Tgi) 


]) 


k-1 


+  41018.88  Ce- 


W  (  X.-X  ) 


o  (H 


i“Xo)  ,k,(P2  expression jk  . 
-H . )  kTrl  14.7  -  1 • 


( .  ■ 


where  P2  expression  is 


P2  =  14.7  + 


0.4  (c  u; 


0.14 


(1  -  C )  (P  -p) 


The  values  of  the  parameters  are  specified  as  follows: 
r  =  the  annual  fixed  charge  factor  =  0.5 
W  =  pounds  of  dry  solids  to  be  treated  per  hour 
=  inlet  solid  moisure  content  =  10% 

=  =  0.111  lb.  water/ lb.  dry  solid 

U  •  v  0 

X  =  outlet  moisure  content  =  1% 
o 

=  q-q  =  -.0101  lb.  water/lb.  dry  solid 
p  =  air  density  =  0.074  lb/cu.  ft. 

ATm  =  (Tpj  ~  Tw)  ~  (TPo  ~  TW) 

TPi  -  Tw 

Lll  m  rp 

APo  -  w 

T  =  wet  -  bulb  air  temperature  =  92°F. 
w 

Tpi  =  inlet  dry  -  bulb  temperature  -  200  F, 

TDq  =  outlet  dry  -  bulb  temperature  is  obtained  according 
to  the  values  assumed  for  Hq  (moisture  content 
of  leaving  air)  by  using  air  -  water  humidity 
charts. 

Table  1  contains  the  different  assumed  values  that 
will  be  used  for  running  the  model. 


•  •  O ' 


Table  1 


Values  of  Outlet  Moisture  Content  of  Air 


Relative 

humidity 

of 

exit  air 

Ho 

.  lb.  water  . 

1  lb. dry  air 

TDo 

(°FJ 

A  Tm 
(°FJ 

90 

0.0325 

95 

29.3 

60 

0.0300 

105  * 

44.8 

40 

0.0275 

115 

55.0 

20 

0 

0.0230 

132 

69.0 

10 

0.0200 

148 

80.0 

5 

0.0135 

174 

93.6 

=  outlet  air  temperature  =  200  F. 

=  inlet  air.  temperature  =  70°F. 

=  moisture  content  of  entering  air  =  0.0008  lb.  water 

lb.  dry  gas 

=  the  overall  heat  -  transfer  coefficient  for 

the  heater  =  5  Bs.tuy'(brJ  (sq.ftj  (°F.) 

=  jVVl  -  (VT.,o).  (250-70)  -  (250-200),  102oF. 


Ts  -  Tqi 


Ln  T  -  T 

s  go 


250-70 

250-200 


(iSzi) 

1  k  ' 


steam  temperature  =  250  F. 
specific  heat  ratio  for  air  =  1.4 
-  2  9 

l  1.4  ' 

solid  density  =  100  lb/cu.ft. 

bed  voidage  =  0.40 

cost  of  steam  $/10^  B.  t.u.  =  2.4 
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It  is  clear  that  some  cost  terms  are  a  direct  function 
of  the  gas  velocity  u  as  in  case  of  dryer  capital  cost, 
compressor  capital  cost,  and  electrical  power  consumption 
cost  while  other  cost  terms  are  independent  upon  gas 
velocity  as  the  heater  capital  cost  and  steam  cost.  For 
stipulated  drying  process  conditions,  the  economical 
gas  velocity  is  obtained  by: 

Minimizing:  Total  cost  =  G(u) 


subject  to 


“mfi  u 


HP 


where  u  ^  =  the  minimum  fluidization  velocity 
uf  =  the  terminal  velocity 
The  optimization  technique  used  in  this  work  is  a 
region  elimination  method  which  narrow  successively  the 
region  contain  the  optimal  solution.  Figure  6  shows  the 

relative  efficiency  of  various  search  optimization  methods. 

\ 

Fibonacci  search  method  is  used  in  this  work  which  is  the 
efficient  one.  The  results  will  be  presented  as  an 
optimal  range  rather  than  a  fixed  value  (9).  The 
Fibonacci  algorithm  is  described  in  Appendix  A. 
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FRACTION  OF  ORIGINAL  INTERVAL  WITHIN  WHICH  OPTIMUM 
LIES  AFTER  PERFORMING  N  EXPERIMENTS.  F 


NIFORM  DICHOTOMOUS 


,-2 1  F  •  O.OI 


SEQUENTIAL  DICHOTOMOUS 


5  PT  RUDD  ft  WATSON(20t 


GOLDEN  MEAN  (I) 


4  8  12  IS  20  24 

NUMBER  OF  EXPERIMENTS,  N 


Figure  6  Relative  efficiency  of  various  search  methods, 
Happel  and  Jordan  (7) 
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IV.  PRESENTATION  AND  DISCUSSION 


OF  RESULTS 

The  model  developed  in  the  previous  chapter  is  used 
to  study  a  fluidized  bed  drying  process.  The  optimal  design 
of  a  fluidized  bed  drying  process  is  affected  by  a  number 
of  factors#  most  important  by  the  fluidized  velocity  and 
the  amount  of  solid  to  be  dried.  Each  component  of  the 
cost  model  is  a  function  of  one  or  both  of  these  factors# 
and  we  are  primarily  interested  in  the  sentitivity  of  cost 
to  changes  in  these  factors.  Other  factors  which  have 
a  minor  effect  on  costs  will  also  be  studied.  The  analysis 
will  address  both  the  analytical  and  technical  reasons 
for  the  cost  behavior. 

This  chapter  will  be  divided  into  three  sections. 

The  first  section  will  analyze  the  behavior  of  the  total 
cost.  The  behavior  of  component  costs  will  be  examined 
in  section  two.  And  finally#  the  third  section  will  review 
the  general  finding  discuss  their  validation#  and  offer 
recommendations  for  further  model  development. 

Since  the  minimum  fluidization  process  is  determined 

by  the  specification  of  the  fluidized  bed  drying  operation# 

a  range  of  this  velocity  is  chosen  according  to  practical 

application  and  examples  in  Leva  (11)  or  Kunii  and 

Levenspiel  (8)  .  This  range  is  from  200  (ft/hx;)  to  1000 

(ft»/hr>)  and  the  base  case  used  for  comparisons  will  be  at 

600  (ft/hrj  .  The  solid  flow  rate  for  the  base  case  will 

be  kept  at  1000  (lby'hrj  and  the  moisture  content  will  be 
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at  0.0325  ( lb.  water/ lb. dry  air) 


Other  factors  which  are  held  constant  in  the  base 
case  and  can  also  affect  capital  cost,  power  cost,  or  both 
are  as  follows: 

1.  Density  of  solid  material  =  100  (lb/cu.ft.) 

2.  Moisture  content  of  outlet  air  =  0.0325  (lb. water) 

(lb. dry  air) 

3.  Electricity  cost  =  .24  ($  /Kw.  hr.) 

4.  Steam  cost  =  2.4  ($  /10^  Bjtu) 

5.  Annual  fixed  charge  factor  =  0.5 

The  sensitivity  of  results  to  these  assumptions  will 
be  examined  in  subsequent  sections  of  this  chapter. 

The  relative  importance  of  equipment  and  power 
components  in  the  total  cost  of  the  drying  operation  are 
shown  in  Table  2.  The  cost  values  are  reported  for  three 
different  optimum  fluidized  velocity  levels.  The  total 
cost  is  1.1  million,  for  each  of  the  dryer  designs.  The 
following  points  emerge  with  respect  to  total  system  cost. 

1.  The  most  important  element  in  cost  is  the  compressor 

It  represents  62.2  -  66.3  per  cent  of  the  total 

cost. 

2.  The  second  most  important  component  of  the  cost 
is  the  dryer  which  represents  32.2  -  36.5  per 
cent  of  total  cost. 

3.  Heater  costs  are  independent  of  fluidized  velocity 
and  represent  a  very  small  percentage  (0.087)  of 


total  cost 


4.  The  three  equipment  components  -  compressor,  dryer, 
and  heater  -  account  for  98.7  per  cent  of  the 
total  cost  while  the  power  consumption  cost  is 
1.3  per  cent  of  the  total  cost. 

Table  2 

Percentage  of  Equipment  Capital  Cost 
and  Power  Consumption  Cost 


Operating 
Fluidized 
Velocity 
(ft/ hr.) 

Cost 

Component  Cost 

Value  x  10^($/year) 

Percentage 

from 

Total  Cost 
(%) 

Dryer 

405.1 

36.5 

Heater 

.9612 

0.087 

10,800 

Compressor 

689.6 

62.2 

Steam 

5.879 

0.53 

Electricity 

7.798 

0.70 

Total  Cost 

1,109 

100.0 

Dryer 

379.9 

34.3 

Heater 

.9612 

0.087 

12,000 

Compressor 

713.0 

64.4 

Steam 

5.879 

0.53 

Electricity 

8.062 

0.73 

Total 

1,108 

100.0 

Dryer 

358.6 

32.4 

Heater 

.9612 

0.087 

13,200 

Compressor 

734.7 

66.3 

Steam 

5.879 

0.53 

Electricity 

8.308 

0.75 

Total 

1,108 

100.00 
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5.  For  a  given  dryer  efficiency,  higher  operating 
velocities  call  for  a  substitution  of  compressor 
horsepower  for  dryer  peripheral  area.  As  shown 

in  Table  2  an  increase  in  the  operating  fluidized 
velocity  from  10,800  (ft/hrj  to  13,200  (ft/hrj 
raises  the  compressor  cost  while  the  dryer  cost 
falls.  The  former  occurs  because  the  compressor 
cost  is  a  direct  function  of  its  horsepower 
which  increases  with  an  increase  in  the  pressure 
required  to  increase  the  fluidization  velocity. 

The  latter  occurs  due  to  the  decrease  in  the 
required  dryer  peripheral  area.  The  decrease  in 
peripheral  area  results  from  the  increase  of  the 
air  flow  rate  required  to  remove  a  fixed  moisture 
content  from  a  given  constant  solid  material. 

This  means  that  for  a  given  drying  operation,  an 
increase  in  fluidization  velocity  will  increase 
the  dry  efficiency  if  we  hold  the  dryer  peripheral 
area  constant  and  consequently  increase  the 
compressor  cost.  However,  if  we  want  the  same 
drying  efficiency  with  velocity  increase  a 
decrease  in  drying  peripheral  area  will  be  feasible 
which  means  a  decrease  in  the  dryer  cost. 

6.  Steam  cost,  which  represents  about  42  per  cent  of 
the  total  power  consumption  cost  is  independent 
of  fluidized  velocity. 
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7.  Electrical  power  consumption#  varies  directly  with 
velocity  and  represent  58  per  cent  of  the  total 


power  cost. 

A.  Capital  and  Power  Cost 

Table  3  and  Figure  7  illustrate  the  sensitivity  of 
total  cost  to  variation  in  velocity.  Two  points  emerge  in 
the  analysis  of  cost  and  velocity.  First#  the  optimum 
velocity  is  independent  of  the  minimum  fluidization  velocity. 
As  shown  in  Table  3  minimum  fluidization  velocity  is 
changed  from  200  (ft. /hr.)  to  1000  (ft. /hr.)  with  the 
results  that  the  optimal  velocity  is  practically  the  same 
and  is  near  the  value  12#000  (ft. /hr.).  Since  the  optimal 
velocity  is  the  same  the  optimal  equipment  design  is  the 
same  and  minimum  system  cost  is  given  at  about  $1100  thousand 
as  shown  in  Table  3.  Figure  7  shows  the  case  of  minimum 
fluidization  velocity  =  600  (ft. /hr.) .  x 

The  results  of  Table  3  will  hold  so  long  as  the 
operating  fluidization  range#  which  is  bounded  by  the  minimum 
fluidization  velocity  and  the  terminal  fluidization  velocity# 
includes  the  optimal  velocity  value.  When  this  is  not 
the  case#  the  minimum  cost  will  occur  at  the  terminal  velocity 
This  fact  is  shown  in  Figures  8  and  9  for  minimum  fluidization 
velocity  values  of  30  and  50  (ft. /hr.)  In  the  cases  optimal 
velocity  is  outside  the  operating  velocity  constant  ratio 
of  the  model.  Figure  7  shows  how  cost  behave  for  the  base 
case  with  minimum  fluidization  velocity  of  600  (ft. /hr.). 


The  optimal  velocity  is  12/202  (ft. /hr.).  Cost  rise 
sharply  for  velocity  values  below  the  optimum  and  rises 
gradually  as  velocity  is  increased  by  and  the  optimum 
range . 

The  extent  of  the  cost  penalty  for  operation  outside 
the  optimum  velocity  range  is  illustrated  in  Table  4. 

Table  4  shows  departure  from  the  optimum  velocity, 

¥ 

12,202  (ft. /hr. ),  and  the  cost  penalties  which  these 
departures  impose  on  the  system.  For  example,  if  the 
system  were  designed  to  operate  at  minimum  fluidization 
velocity  if  600  ft. /hr.,  system  costs  would  be  $2.6  million  or 
137.30  per  cent  more  than  would  be  incurred  if  the  system 
subject  to  this  minimum  velocity  were  designed  to  operate 
at  the  optimum  fluidization  velocity  of  12,202.  If  the 
system  were  designed  to  one-half  the  optimum  velocity 

(i.e.  velocity  of  6,100  (ft. /hr.),  cost  would  be  $1.16  million 

\ 

or  4.7  per  cent  greater  than  the  minimum  system  cost. 

A  similar  cost  penalty  occurs  with  systems  designed 
to  operate  at  fluidization  velocities  in  excess  of  the 
optimum.  For  example  at  twice  the  velocity  level,  the 
system  design  would  cost  $1.15  million  or  4.15  per  cent 
more  than  the  optimal  design.  At  very  high  velocities, 
the  cost  penalty  is  not  nearly  as  great  as  it  is  for 
very  low  fluidization  velocities.  For  example,  at  the 
terminal  velocity  of  54,600  (ft. /hr.)  the  system  design 
would  cost  202  thousands  ($/yr<)  or  about  18.23  per  cent  more 
than  the  minimum  cost  design. 


The  greater  cost  penalty  of  suboptimal  velocity 
designs  reflects  the  fact  that  the  compressor  cost  is 
minimized  at  the  minimum  fluidization  velocity  while  the 
dryer  cost  is  minimized  at  the  terminal  velocity. 

Table  3 

Capital  and  Power  Cost 
at  Different  Minimum  Fluidization  Velooities 


Minimum  Optimal  Veloc 

Fluidization  Range  (ft/hrj 

Velocity  _ _ 

(fb/hr.) 


Optimal  Velocity  Minimum  Cost  Range 

Range  (ft/hrj  x  10^  ($/yr;) 

o  From  To  ” 


12,231 

12,184 

1,107.77 

1,107.77 

12,202 

12,345 

1,107.77 

1,107.78 

12,220 

11,981 

1,107.76 

1,107.80 

FIGURE <7).  CAPITAL  AND  POWER  COST 


3 


Table  4 

Cost  Penalties  for  Operating 
Outside  the  Optimal  Velocity  Range  (12,202  ftVhn) 


Case 

Fluidized 

Velocity 

(ft>/hrO 

Difference 
from  Optimal 
Velocity 
(f  fc/hrj 

Cost  .103 
($/yzJ 

Difference 

from 

Minimum 

cost 

x  103  ($/hc) 

Per  cent 
Increase 

1 

600 

-11602 

2,629 

1,521 

137.30 

2 

6100 

-6102 

1,160 

52 

4.7 

3 

13800 

1598 

1,109 

1 

0.09 

4 

20400 

8198 

1,133 

25 

2.26 

5 

24400 

.12198 

1,154 

46 

4.15 

6 

33600 

21398 

1,203 

95 

8.57 

7 

40200 

27998 

1,239 

131 

11.82 

8 

46800 

34598 

1,272 

164 

14.80 

9 

54600 

42398 

1,310 

202 

18.23 
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2.  Solid  Flow  Rate 


A  change  in  the  solid  flow  rate  will  affect  each 
component  of  the  cast  model.  The  reasons  for  this  ares 

1.  Dryer  costs  will  increase  with  the  solid  flow 
rate  increase  due  to  the  change  in  its  design. 

The  dryer  needed  for  higher  solid  flow  rate  will 
have  a  greater  volume  and  a  greater  diameter, 
consequently  this  greater  diameter  increases  the 
peripheral  area  which  is  the  determining  factor 
in  dryer  capital  costs. 

2.  Compressor  costs  will  also  increase  since  a 
higher  pressure  drop  will  be  needed  to  fluidize 
the  solid  material.  This  requires  a  higher  horse¬ 
power  which  leads  to  a  higher  capital  cost. 

3.  The  heater  cost  will  increase  with  the  solid 

flow  rate  increase  due  to  the  increase  in  the 

\ 

total  moisture  of  the  solid  material.  This 
requires  more  heating  load  which  increases  the 
peripheral  area  of  the  heater  and  consequently 
its  capital  cost. 

4.  Power  consumption  increases  with  a  greater  solid 
flow  rate  due  to  the  increase  in  pressure  drop 
needed  for  the  fluidization  process.  This 
increases  the  electrical  power  consumption. 

Also  the  increased  heating  load  required  to  remove 
the  increased  amount  of  solid  moisture  content 

increase  the  steam  consumption. 
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Table  5  shows  the  optimal  velocity  change  and  the 
minimum  cost  change  due  to  the  changes  in  the  solid  flow 
rate.  If  the  system  were  designed  to  operate  at  twice 
the  solid  flow  rate,  a  decrease  in  optimal  fluidization 
velocity  by  about  14  per  cent  which  causes  an  increase 
in  the  minimum  cost  of  1  million  or  90.8  per  cent. 

At  high  flow  rates  of  7000  (lbVhrJ  the  optimal  velocity 
decreases  to  about  one-half  the  optimum  velocity  of  the 
base  case  (12,202  (ft-/h&)  while  the  cost  increases  by 
about  5.7  million  or  513  per  cent. 
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FIGURE <10>.  CAPITAL  AND  POWER  COST 


3.  Other  Auxiliary  Factors 

Other  factors  that  have  a  slight  influence  on  the  cost 
model  can  be  divided  as  major  factors  and  minor  factors. 

Major  factors  are  the  density  of  solid  material  and  the 
moisture  content  of  the  outlet  air  while  minor  factors  are 
the  steam  cost/  the  electricity  cost,  and  the  annual  fixed- 
change  factor.  Effects  of  variation  in  these  factors  on 
model  results  are  shown  in  Table  5.  The  essential  effects 
are  as  follows. 

1.  Variation  in  solid  density  will  require  adjustments 
in  the  pressure  drop  and  consequently  in  the  com¬ 
pressor  and  its  capital  cost  and  the  electrical 
power  consumption  of  the  process.  For  example,  an 
increase  in  solid  density  to  150  (lb/cu. ffaj  results 
in  a  decrease  in  the  optimum  velocity  range,  a 
system  design  which  includes  a  compressor  of  greater 
horsepower  and  an  increase  in  system  cost  of  about 
$300  thousands. 

2.  Variation  in  moisture  content  of  outlet  air  from 
0.0325  to  0.023  { lb. water/ lb, dry  air)  has  a  great 
effect  on  the  optimal  fluidization  velocity  and  a 
decrease  in  system  cost  of  about  $500  thousands. 

The  optimal  fluidization  velocity  reaches  the 
terminal  velocity  in  this  case.  Decrease  the  outlet 
moisture  content  means  increasing  in  dry-bulb 
temperature  of  air  as  indicated  in  air-water 


humidity  charts  (14).  The  effect  of  decreasing 
outlet  moisture  content  of  air  is  an  increase  in 
the  optimal  fluidization  velocity  and  consequently 
an  increase  in  the  compressor  cost.  On  the  other 
hand,  the  dry-bulb  temperature  is  increased  which 
increases  the  logarithmic  mean  temperature  and 
consequently  the  dryer  peripheral  area  is  decreased 
which  decreases  the  dryer  cost.  This  decrease  in 
dryer  cost  is  higher  than  the  increase  in  compressor 
cost,  so  the  net  result  of  this  case  is  a  decrease 
in  total  cost. 

Variation  in  the  steam  cost  from  4.8  to  1.6 
($/10^  BJbuJ  has  no  effect  upon  the  optimal 
velocity  range.  However  it  increases  or  decreases 
slightly  the  capital  and  power  cost. 

Variation  in  the  electricity  cost  from  0.48  to 
0.16  ($/Kw. hrj  has  no  effect  upon  the  optimal 
velocity  range.  The  capital  and  power  cost  will 
increase  or  decrease  slightly  more  than  was  the 
case  of  variation  in  the  steam  cost. 

Variation  in  the  annual  fixed  change  factor  from 
0.4  to  0.6  will  have  no  impact  on  the  optimal 
velocity  range.  Of  course  variation  in  this  change 
will  increase  or  decrease  total  cost. 


Table  6 

Capital  and  Power  Cost  at  Different 
Auxiliary  Factors 


Factor  to  be  Changed 


Optimal  Velocity 
Range  (ft. /hr.) 


Minimum  Cost 
Range  x103  ($/yrJ 


Value 


Base  Case 

12,202 

12,345 

1,107.77 

1,107.78 

Solid 
density 
(lb/cu.  f  tj 

150 

8,541 

6,953 

1,405.99 

1,405.65 

Moisture 
content 
of  the 
outlet 
air 

,1b. water/  . 
'lb. dry  air' 

0.023 

54,457 

54,600 

\ 

1,056.14 

1,054.96 

Steam  Cost 
($/106  B.tu.) 

Cs 

1.6 

12,202 

12,345 

1,105.81 

1,105.82 

Electricity 

Cost 

($/Kw.hrJ 

ce 

.16 

12,202 

12,345 

1,105.07 

1,105.07 

Annual 

Fixed 

Change 

Factor  R 

0.6 

12,202 

12,345 

1,326.53 

1,326.53 

Separate  Cost  Items 


In  this  section  the  behavior  of  each  of  the  component 


cost  items  is  briefly  reviewed. 


li  Capital  Cost 

Capital  cost  represents  98.7  per  cent  of  the  total 
cost  and  it  is  obvious  that  the  behavior  of  the  total 
cost  is  mainly  due  to  the  behavior  of  capital  cost.  Table  7 
and  Figure  11  illustrate  the  behavior  of  capital  cost.  A 
comparison  of  Table  7  with  Table  3  of  the  previous  section 
shows  that  the  optimal  velocity  range  is  the  same,  and  the 
slight  difference  in  the  minimum  cost  range  is  due  to  the 
presence  of  power  cost  in  Table  3. 

Capital  cost  components  are  the  dryer  cost,  the 
compressor  cost,  and  the  heat  exchanger  cost.  Examination  of 
the  component  capital  cost  shows  the  following: 

1.  As  shown  in  Table  8  and  Figure  12  the  minimum  cost 
dryer  design  is  at  the  upper  limit  of  operating 
velocity  which  is  the  terminal  velocity. 

2.  Compressor  capital  cost  increases  as  fluidized 
velocity  increases  as  can  be  seen  in  Table  9  and 
Figure  13. 

3.  The  heat  exchanger  capital  cost  is  independent  of 
velocity,  has  a  fixed  value  at  961  ($/yrO ,  and  is 
very  low  compared  with  compressor  or  dryer  capital 
costs  representing  about  0.1  per  cent  of  total 
capital  cost. 


The  optimal  drying  operation  design  results  in  a  kind 
of  compromise  between  the  dryer  and  compressor  costs 
corresponding  to  an  optimal  velocity  intermediate  between  the 
minimum  and  terminal  fluidization  velocities  of  the  system. 

Table  7 
Capital  Cost 

and  Different  Minimum  Fluidization  Velocities 


Case 

Minimum 

Fluidization 

Optimal  Velocity 
Range  (ft>/hr0 

Minimum  Cost  Range 
Xl03  ]($/yrJ 

velocity  uty  nr.; 

From 

To 

From 

To 

1 

200 

12,375 

12,327 

1,093.77 

1,093.77 

2 

600 

12,345 

12,489 

1,093.77 

1,093.78 

3 

1000 

12,459 

12,220 

1,093.77 

1,093.78 

Table  8 
Dryer  Cost 

At  Different  Minimum  Fluidization  Velocities 


Case 

Minimum 
Fluidization 
Velocity  (ft/hrj__ 

Optimal  Velocity 
Range  (ft/hrj 

Mimimum  Cost  Range 
Xl03  ($/yri) 

From 

To 

From 

To 

1 

200 

18,152 

18,200 

295.42 

294.94 

2 

600 

54,457 

54,600 

151.48 

151.23 

3 

1000 

90,761 

91,000 

111.09 

110.86 

Table  9 


Compressor  Cost 

at  Different  Minimum  Fluidization  Velocities 


Case 

Minimum 
Fluidization 
Velocity  (ftVhrj 

Optimal  Velocity 
Range  (fb/hrj 

Minimum  Cost  Range 
yl03  { $  /  y  rj 

From 

TO 

From 

TO 

1 

200 

248 

200 

201.71 

187.84 

2 

600 

743 

600 

290.18 

270.41 

3 

1000 

1239 

1000 

343.19 

319.92 

FIGURE <1 1). CAPITAL  COST 


VELOCITY  ^ 

FIGURE (12).  ORYER  COST 


y 

2.  Power  Consumption  Cost 

The  power  consumption  cost  represents  1.3  per  cent  of 
the  total  cost.  The  components  of  power  cost  are  the 
steam  cost  and  the  electricity  cost.  The  behavior  of  power 
consumption  cost  and  its  components  is  shown  in  Figures 
15/  16  and  17.  As  can  be  seen: 

1.  The  power  consumption  cost  increases  with  velocity 

2.  Steam  cost  is  independent  of  velocity,  and 

3.  The  electrical  consumption  cost  is  responsible  for 
the  direct  relationship  of  power  cost  to  velocity. 

The  sensitivity  of  power  cost  to  changing  values  of 
steam  C  and  electricity  C  is  shown  in  Table  10.  As  can 

S  G 

be  seen  a  doubling  or  halving  of  evergy  price  raises 
(reduces)  power  costs  30  (30)  per  cent.  Optimal  velocity 
and  dryer  design  is  little  affected  by  small  variation  as 
noted  above  in  Section  A. 
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Table  10 

Power  Consumption  Cost 
at  Different  Steam  and  Electricity  Prices 


Ce 

Cs 

Optimal  Velocity 
Range  (fb/hrj 

Minimum  Cost  Range 
.103  ($/yiJ 

From 

TO 

From 

To 

.16 

1.6 

743 

600 

\ 

6.107 

5.958 

.16 

2.4 

743 

600 

8.067 

7.918 

.24 

2.4 

743 

600 

9.160 

8.934 

.24 

1.6 

743 

600 

7.201 

6.978 

FIGURE (17). STEAM  CONSUMPTION  COST 


C.  Summary,  Conclusions  and  Recommendations 

The  cost  model  developed  and  analyzed  in  this  study 
can  be  employed  to  identify  optimal  fluidization  velocities 
and  their  associated  minimum  cost  equipment  designs.  The 
model  results  can  be  summarized  as  follows: 

1.  Capital  costs  are  very  high  compared  to  power  costs. 
The  former  represents  98.7  per  cent  of  the  total 
cost  of  an  optimum  dryer  design  configuration. 

2.  Dryer  capital  cost  varies  inversely  with  fluid 
velocity  and  represents  33  -  37  per  cent  of  the 
capital  cost  of  an  optimal  design  configuration. 

3.  Compressor  costs  vary  directly  with  fluid  velocity 
and  represent  63  -  67  per  cent  of  the  capital  cost 
of  an  optimal  design. 

4.  The  optimal  capital  design  requires  a  dryer,  compressor 
and  heater  configuration  which  has  an  optimal 
fluidized  velocity  which  is  intermediate  to  the 
minimum  fluidization  velocity  (which  gives  lowest 
compressor  cost)  and  the  terminal  velocity  (which 
gives  lowest  dryer  costs) . 

5.  Heater  requirements  increases  the  capital  costs 
of  the  optimum  design  at  the  indicated  optimal 
velocity  level. 

6.  The  power  cost  is  relatively  unimportant  representing 


1.3  per  cent  of  total  design  cost.  Of  this,  steam 
requirements  represents  42  per  cent  of  the  power 


cost,  while  electricity  requirements  represents 
58  per  cent  of  total  power  cost.  Only  the  latter 
varies  with  fluidization  velocity,  increasing  as 
the  pressure  drop  increases. 

7.  Equipment  designs  outside  the  optimal  fluidization 
velocity  range  impose  cost  penalties.  Cost 
penalties  can  be  held  to  less  than  4  per  cent  of 
system  costs  so  long  as  design  configurations 
correspond  to  fluid  velocities  no  less  than  one  half 
the  optimal  velocity  nor  more  than  double  the 
optimal  velocity.  Very  severe  cost  penalties  are 
imposed  on  system  designed  at  or  close  to  minimum 
fluidization  velocity  levels. 

A  brief  check  in  the  validity  of  the  design  model  is 
possible  by  checking  two  key  aspects  of  the  model  against 

reported  real  world  drying  operations.  First  the  real  world 

\ 

examples  of  fluidized  velocity  drying  operation  report 
velocity  levels  which  lie  inside  the  optimal  velocity 
range  identified  by  the  model  (11) ,  (8) . 

Second  the  height/diameter  ratio  reported  to  be  most 
efficient  in  commercial  practice  lies  between  4  and  10 
(15)  #  Solutions  of  the  model  for  a  solid  flow  rate  from 
7,800  to  10,000  (lb/hrj  and  fluidized  velocity  in  the 
optimal  range  5,000  to  8,000  (ft/hrj  yields  a  ratio  of 
height/diameter  of  4  to  5.  These  limited  information  on 
operating  systems  is  agreed  with  the  results  of  our  model. 
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The  model  developed  proposed  in  this  study  represents 
a  first  step  in  the  development  of  a  mathematical  model 
for  use  in  the  design  and  study  of  f luidized-bed  drying 
operations.  Model  modification  to  include  the  height- 
diameter  ratio  as  a  constraint  in  the  design  process  would 
be  a  useful  extension  of  this  model.  Also  the  mechanisms 
of  f luidized-bed  bubble  formation  and/or  channeling  and 
the  mechanism  effect  on  heat  transfer  might  be  usefully 
considered  in  future  development  of  the  model. 
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APPENDIX  A 


o* 


Fibonacci  Algorithm 

This  algorithm  finds  the  minimum  of  a  single  variable, 
nonlinear  function  subject  to  constraints. 

Minimize  F  (X) 

subject  to  a.^  <  x  <  b^ 

The  upper  and  lower  bounds,  bj^  and  a^^  are  constants. 

The  procedure  is  an  interval  elimination  search 
method.  Thus,  starting  with  the  original  boundaries  on  the 
independent  variable,  the  interval  in  which  the  optimum  value 
of  the  function  occurs  is  reduced  to  some  final  value, 
the  magnitude  of  which  depends  on  the  desired  accuracy. 

The  location  of  points  for  functions  is  based  on  the  use  of 
positive  integers  known  as  the  Fibonacci  numbers.  No 
derivatives  are  required.  A  specification  of  the  desired 
accuracy  will  determine  the  number  of  function  evaluations. 

A  unimodal  function  is  assumed.  The  algorithm  proceeds 
as  follows: 

1.  Designate  the  original  search  interval  as  with 
boundaries  a^  and  b^. 

2.  Predetermine  the  desired  accuracy  a  and  thus  the 


number,  N,  of  required  Fibonacci  numbers. 


F 

n 
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where  F  is  called  a  Fibonacci  number 
n 

3.  Place  the  first  two  points,  X1  and  X2  (X^^  <  X2> 

within  at  a  distance  1^  from  each  boundary. 

F 

r  -  N-2 

*1  F  L1 

A  N  A 

Xi  =  ai  + 

x2  =  bl  ‘  ll 

4.  Evaluate  the  objective  function  at  X^  and  X2» 
Designate  the  functions  as  FfX^  and  F(X2>. 
Narrow  the  search  interval  as  follows: 


ai  <  X*  <  x2 


for  F(x1)  <  f(x2) 


Xx  <  X*  <  bx 


for  FfX^  >  F(X2) 


Where  X*  is  the  location  of  the  optimum.  The  new  search 


interval  is  given  by 

L  -  Fn~1 
"2  F 

*  *n 


•  li  =  li  - 


with  boundaries  a2  and  b2 

5.  Place  the  third  point  in  the  new  L2  subinterval, 

symmetric  about  the  remaining  point, 

F 

i  =  A"3  l 

1  *N-1 

X3  —  ^2  +  12  or  b2  —  12 

6.  Evaluate  the  objective  function  F(X.j)  ,  compare  with 
the  function  for  the  point  remaining  in  the 


interval  and  reduce  the  interval  to 
F 

t  _  N-2  _ 

L2  "  Fn  L1  =  L2  "  h 
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7.  The  process  is  continued  per  the  preceeding  rules  for 


N  evaluations, 
v  =  FN  ~  (k+1) 

k  FN  -  (k-1) 

XL+1  =  ak  +  Lk  or 

L.  =.  V*.  ■tfc.-.M. 


‘k 

L 


(symmetric  about  mid  point) 
k-1  “  1k-l 


After  N-l  evaluations  and  discarding  the  appropriate 
interval  at  each  step,  the  remaining  point  will  be  precisely 
in  the  center  of  the  remaining  interval.  Thus  point  N  is 
also  at  the  midpoint  and  is  replaced  by  a  point  perturbed 
some  small  distance  e  to  one  side  or  the  other  of  the 
midpoint.  The  objective  function  is  then  evaluated  and  the 
final  interval  where  the  optimum  is  located  is  thus  determined 


APPENDIX  B 


ij 


100=*************************************************************** 

110=*  THIS  A  FORTRAN  CODE  OF  FIBONACCI  SEARCH  ALGORITM  USED  IN  * 

120=*  FINDING  THE  OPTIMUM  FLUIDIZATION  VELOCITY  OF  FLUIDIZED -BED  * 

130=*  DRYINC  OPERATION.  * 

X 40s******************** ********************************  *********** 

150=  PROGRAM  FIBONA 

160=  DIMENSION  FIB(50) 

170=  ALPHA-. 01 

180=  A=6  00 

190=  A1»A 

200=345  B=91 *A 

210=  DEL=B-A 

220=  WRITE ( * , 001 ) 

230=001  FORMAT (//////) 

240=567  FORMAT ( 1 5X , '  *  MINIMUM  FLUIDIZATION  VE LOC ITY= ' , F5 . 0 , ’ ( FT/ HR ) * ' ) 

250  =  223  F0RMAT(  22X  ,  '  *  110=  .  03  2  5  ,  TD0  =  9  5  .  ,  W  =  1 000  *') 

260  =  789  FORM AT( 2 2X ,  ' *  CAPITAL  AND  POWER  COST  *') 

270=  WRITE ( *  ,  567  )A 

280=  WRITE ( * , 223 ) 

290=  WRITE ( * , 789 ) 

300-C 

3 1 0=C  DEFINE  THE  FIRST  THREE  FIBONACCI  NUMBERS 

320-C 

330=  FIB0-1.0 

340=  FIB(1)  =  1 .0 

350=  FIB(2)=2.0 

360-C 

370-C  CALCULATE  THE  REMAINING  FIBONACCI  NUMBERS 

380-C 

390=5  BB=1. O/ALPHA 

400=  IF(BB-2. 0)10, 10,11 

410=10  GO  TO  14 

420*11  CONTINUE 

430=  J J=2 

440=12  JJ-JJ+1 

450=  FIB( JJ)=FIB( JJ-1 )  +  FIB( JJ-2  ) 

460=  CC  =  F IB ( J  J ) 

470=  IF(CC-BB) 13, 15, 15 

480=13  GO  TO  12 

490=14  WRITE ( * , 002 ) 

500=002  F0RMAT(/// , 10X, 'ACCURACY  SPECIFIED  IN  FUNC  NOT  SUFFICIENT', 
510=  1//.10X, 'PROGRAM  RESET  ALPH A , ALPHA- . 01 ' ) 

520=  ALPHA-. 01 


GO  TO  5 


530- 
540-C 

550-C  FIRST  STEP  IN  THE  TABLEAU 

560-C 

570-15  1-0 

580-  KK-JJ-2 

590-  IK-JJ-2 

600-  BL-B-A 

610-  ALL-FIB(IK)*BL/FIB( JJ) 

620-  W-A+ALL 

630-  V-B-ALL 

640-  CALL  FUNC(W,T) 

650-  CALL  FUNC (  V  ,  U  ) 

660-  JK-1 

670-  WRITE ( * , 555 ) 

680-  WRITE (*,777) 

690-  WRITE ( *  ,003  ) 

700-003  FORMAT (  ,**,,,K•,,  * » , 4X , ' LK • , 4X , • * • , 4X , • AK • , 4X , ' * ' , 

710-  +4X,  '3K  '  ,4X,  ' , 3X, »LLK* ,4X, 1 '  VELOCITY  TOTAL  COST 

72  0-777  FORMAT (IX,  • * '  , 2X , *  * '  , 1 OX , 1  * '  , I  OX , ’ * • , IOX ,  ' *  •  , 10X , • * ' , 

730=  +  IOX, •*» , 1 2X , '*' ) 

740-  WRITE (*,777) 

750-  WRITE ( *,555) 

760-  WRITE (*,777) 

770-  WRITE (*, 004 )JK,BL, A, B, ALL,W,T 

780-  WRITE ( *,777) 

790-  WF ITE ( *,555) 

800-  WRITE ( *,777) 

810-004  FORMAT (  '**,,I2,,*,,E10.4,,*',E10.4,,*,,E10.4,'*,,E10.4, 

820-  I  ,E10.4, • *' ,E12.6, **'  ) 

830-006  FORM AT (5 IX, ,*',E10.4,'*,,E10.4) 

840-C 

850-C  SUCCEEDING  STEPS  IN  THE  TABLEAU 

860-C 

870-  IK-IK-1  • 

880-  JJ-JJ-1 

890-  DO  70  I-l.KK 

900-  IF(U-T)20, 20, 22 

910-20  A-A+ALL 

920-  BL-B-A 

930-  W-V 

940-  CALL  FUNC ( W , T ) 

950-  ALL-FIB(IK)*BL/FIB(JJ) 

960-  V-B-ALL 

970-  CALL  FUNC( V , U) 

980-  II-I+l 

990-  IK-IK-1 

1000-  JJ-JJ-1 

1010-  IF(IK-1)28, 29, 29 

1020-28  IK-1 
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1030-29  CONTINUE 

1040-  WRITE(*,004)II ,BL,A,B, ALL,W,T 

1050-  WRITE (*,777) 

1060-  WRITE ( *,555) 

1070-  WRITE ( *,777) 

1080-  GO  TO  70 

1090-22  B-B-ALL 

1100-  BL-B-A 

1110-  V-W 

1120-  CALL  FUNC (  V  ,  U  ) 

1130-  ALL»FIB(IK)*BL/FIB(JJ) 

1140-  W-A+ALL 

1150-  CALL  FUNC( W,T) 

1160-  II-I+l 

1170-  IK-IK-1 

1180-  JJ-JJ-1 

1190-  IF(IK-1)30,31,31 

1200-30  IK-1 

1210-31  CONTINUE 

1220-  WRI TE(*, 004)11, BL, A, B, ALL ,  V ,  U 

1230-  WRITE( * , 555  ) 

1240-555  F0RMAT(1X,72( ’*' )) 

1250-  GO  TO  70 

1260-70  CONTINUE 

1270-C 

1280-C  CALCULATION  OF  THE  FINAL  RANGE  OF  THE  DEPENDENT  VARIABLE 

1290-C 

1300-  EPS- . 001 *W 

1310-  DL-EPS+W 

1320-  CALL  FUNC(DL.YL) 

1330-  IF(YL-T)80,80,81 

1340-80  CALL  FUNC(B.BF) 

1350-  WRITE (*,333) 

1360-333  F  ORM AT (////, 2X) 

1370-  WRITE( * , 567 )Al 

1380-  WRITE ( *,223) 

1390-  WRITE ( * , 789 ) 

1400-  WRITE(*,666) 

1410-  WRITE ( * , 444 ) 

1420-4  44  FORMAT (IX, ' * ' , 24X , ' * ' , 1 5X , • * ' , 1 5X , ' *  * ) 

1430-  WRITE(*,007)W,B 

1440-  WRITE C*  , 444 ) 

1450-007  FORMAT( *  **  THE  OPITMAL  VELOCITY  * ' , 2X , E 1 3 . 7 , • * • , 2X , E 1 3 . 7 
1460-  WRITE ( * , 666 ) 

1470-  WRITE ( *  ,  444 ) 

1480-666  F0RMAT(1X,58( *  *  '  )) 

1490-  WRITE ( *,008)T,BF 

1500-  WRITE ( *  ,  444  ) 

1510-008  FORMAT (  '  **  CAPITAL  AND  POWER  COST  *  '  ,  2X ,  F.  1  3 . 7  ,  *  *  '  ,  2X ,  El 3 . 7 
1520-  WRITE ( * , 666  ) 

1530-  WRITE(*, 333) 
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1540-  CO  TO  87 

1550-81  CALL  FUNC(A.AF) 

1560-  WRITE(*,007)W,A 

1570-009  F0RMAT(///,f  THE  FINAL  FEASIBLE  REGION 2X X- ’  , 

1580-  1  E13.7.2P, 'X-',E13.7) 

1590-  WRITE(*,008)T,AF 

1600-  WRITE( * , 333  ) 

1610-017  FORMAT ( / '  WITH  FUNCTION  VALUE S ' , 7X , ’ Y- '  , El 3 . 7 , 2X , ' Y- • , 

1620-  +  E10.4) 

1630-87  ACC- ( W-A) /(DEL) 

1640-999  CONTINUE 

1650-123  END 
1660-C 

1670-  SUBROUTINE  FUNC(X,Y) 

1680-  W-1000.0 

1690-  XI-.lll 

1700-  X0-.0101 

1710-  TDI-200 

1720-  TD0-95 

1730-  TW-92 

1740-  HI-. 0008 

1750-  ROH- .074 

1760-  DELM-( (TDI-TW)-(TDO-TW) ) / ( ALOG(TDI-TW) /(TDO-TW) ) 

1770-  EP-.4 

1780-  ROS-150 

1790-  P2-14.7+.4*( (ROH*X)** . 34 ) * ( 1 -EP ) * ( ROS-ROH ) /DELM 

1800-  HO- .0325 

1810-  UH-5 

1820-  TGO-200 

1830-  TGI -7  0 

1840-  TS-250 

1850-  DELH*((TS-TGI)-(TS-TGO))/(ALOG(TS-TGI)/(TS-TGO)) 

1860-  SPCR-.29 

1870-  R-.5 

1880-  CE-0.24 

1890-  CS-2.4/.1E+07 

1900-  DRCT-  R*(585*((200*W*(XI-XO))/(ROH*X)**.76*DELM)**.8) 

1910-  HTCT-R*(346*(W*(XI-XO)*( .24* (TGO-TGI )+ 

1920-  +  .5*HI*(TG0-TGI))/((II0-HI)*U1I*DELH))**.62) 

1930-  COCT-R*5873 . 6 3* < W* ( XI-XO ) *S PCR* ( (P2/14. 7)* *SPCR-l)/( (R0H*( 

1940-  +  HO-HI) ) )** .8) 

19  50-  STCT-8000*CS*(W*( XI-XO )*( . 24  * (TGO-TGI )+ 

1960-  +  . 5* ( T  GO -TGI ))/(HO-HI)) 

1970-  POCT-4 1018.88*CE*W*(XI-XO)*SPCR*((P2/14.7)**SPCR-1)/ 

1980-  +  ROH* ( HO-H I ) 

1990-  CAPCT-DRCT+HTCT+  STCT+POCT+COCT 

2000-  Y-CAPCT 

2010-  RETURN 

2020-  END 

2030-*E0R 


i 
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